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The surface of a three-dimensional topological electron system often hosts symmetry-protected 
gapless surface states. With the effect of electron interactions, these surface states can be gapped 
out without symmetry breaking by a surface topological order, in which the anyon excitations carry 
anomalous symmetry fractionalization that cannot be realized in a genuine two-dimensional system. 
We show that for a mirror-symmetry-protected topological crystalline insulator with mirror Chern 
number n = 4, its surface can be gapped out by an anomalous Z 2 topological order, where all anyons 
carry mirror-symmetry fractionalization M 2 — —1. The identification of such anomalous crystalline 
symmetry fractionalization implies that in a two-dimensional Z 2 spin liquid the vison excitation 
cannot carry M 2 = — 1 if the spinon carries M 2 = — 1 or a half-integer spin. 
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The advent of topological insulators (TIs) (THU] and 
topological superconductors (TSCs) U has greatly 
broadened our understanding of topological phases in 
quantum systems. While the concepts of TIs and TSCs 
originate from topological band theory of noninteract¬ 
ing electrons or quasiparticles, recent theoretical break¬ 
throughs EHim have found that interactions can, in prin¬ 
ciple, change fundamental properties of these topological 
phases dramatically, thus creating a new dimension to 
explore. In particular, interactions can drive the gapless 
Dirac fermion surface states of three-dimensional (3D) 
TIs and TSCs into topologically ordered phases that are 
gapped and symmetry preserving. Nonetheless, such a 
surface manifests the topological property of the bulk in a 
subtle but unambiguous way: its anyon excitations have 
anomalous symmetry transformation properties, which 
cannot be realized in any two-dimensional (2D) system 
with the same symmetry. 

Given the profound consequences of interactions in 
TIs and TSCs, the effect of interactions in topologi¬ 
cal phases protected by spatial symmetries of crystalline 
solids, commonly referred to as topological crystalline in¬ 
sulators (TCIs) 111 . is now gaining wide attention. A 
wide array of TCI phases with various crystal symme¬ 
tries have been found in the framework of topological 
band theory mug. One class of TCIs has been pre¬ 
dicted and observed in the IV-VI semiconductors SnTe, 
Pbi_a;Sn x Se and Pbi-^Sn^Te nun]. The topologi¬ 
cal nature of these materials is warranted by a partic¬ 
ular mirror symmetry of the underlying rocksalt crystal, 
and is manifested by the presence of topological surface 
states on mirror-symmetric crystal faces. Remarkably, 
there surface states were found to become gapped un¬ 
der structural distortions that break the mirror symme¬ 
try mm, confirming the mechanism of crystalline pro¬ 
tection unique to TCIs |TTj . 

The study of interacting TCIs has just begun. A recent 


work by Isobe and Fu [20] shows that in the presence of 
interactions, the classification of 3D TCIs protected by 
mirror symmetry (i.e., the SnTe class) reduces from being 
characterized by an integer known as the mirror Chern 
number [2T] (hereafter denoted by n ) to its Zg subgroup. 
This implies that interactions can turn the n = 8 surface 
states, which consists of eight copies of 2D massless Dirac 
fermions with the same chirality, into a completely triv¬ 
ial phase that is gapped, mirror symmetric and without 
intrinsic topological order. It remains an open question 
what interactions can do to TCIs with n / 0 mod 8. In 
this work, we take the first step to study strongly inter¬ 
acting TCI surface states for the case n = 4 mod 8. 

Our main result is that the surface of a 3D TCI 
with mirror Chern number n = 4 mod 8 can become 
a gapped and mirror-symmetric state with Z 2 topolog¬ 
ical order. Remarkably, the mirror symmetry acts on 
this state in an anomalous way that all three types 
of anyons carry fractionalized mirror quantum number 
M 2 = -1 (in this Letter we use M to represent the pro¬ 
jective representation of mirror symmetry M acting on 
an anyon), which cannot be realized in a purely 2D sys¬ 
tem. Furthermore, the anomalous mirror-symmetry frac¬ 
tionalization protects a twofold degeneracy between two 
mirror-symmetry-related edges. Such anomalous mirror- 
symmetry fractionalization cannot be realized in a 2D 
system, including a 2D Z 2 spin liquid state [22] , Hence, 
our finding constrains the possible ways of fractionalizing 
the mirror symmetry in a 2D Z 2 spin liquid I23 U 21 I . Brief 
reviews of 3D TCIs, 2D Z 2 spin liquids, and their edge 
theory are available in the Supplemental Material [251 . 

Noninteracting TCIs. We begin by considering nonin¬ 
teracting TCIs protected by the mirror symmetry x -A 
— x . With the mirror symmetry, the extra U(l) symme¬ 
try in a TCI does not change the classification in 3D 
compared to a mirror-protected topological crystalline 
superconductor. Hence, for convenience, we choose a 
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TCI as our starting point, although the U(l) symme¬ 
try plays no role in this work. As we will explain in 
Sec.[n]of the Supplemental Material, in order to produce 
an anomalous Z 2 surface topological order, the mirror 
operation must be defined as a Z 2 symmetry with the 
property M 2 = 1. In our previous works on spin-orbit 
coupled systems, the mirror operation M' acts on elec¬ 
tron’s spin in addition to its spatial coordinate, which 
leads to M' 2 = —1. Nonetheless, one can redefine the 
mirror operation by combining M' with the U(l) sym¬ 
metry of charge conservation c —> ic , which restores the 
property M 2 = 1. We note that without the U(l) sym¬ 
metry only M satisfying M 2 = +1 protects nontrivial 
topological crystalline superconductors. 

The mirror TCIs are classified by the mirror Chern 
number n, defined for single-particle states on the mirror- 
symmetric plane k x = 0 in the 3D Brillouin zone. The 
states with mirror eigenvalues 1 and —1 form two dif¬ 
ferent subspaces, each of which has a Chern number de¬ 
noted by n + and n_ respectively. This leads to two in¬ 
dependent topological invariants for noninteracting sys¬ 
tems with mirror symmetry: the total Chern number 
riT = n++n_ and the mirror Chern number n = n+—n-. 

The TCI with a nontrivial mirror Chern number n has 
gapless surface states consisting of n copies of massless 
Dirac fermions, described by the following surface Hamil¬ 
tonian 

n 

H s = v ^2 i/j' A (k x a y - k y a x )xp A , (1) 

A =1 

where the two-dimensional fermion fields ipA{x,y) trans¬ 
form as the following under mirror operation: 

M : ipA{x, y) -> a x ip A {~x,y). (2) 

The presence of mirror symmetry Q forbids any Dirac 
mass term As a result, the surface states de¬ 

scribed by Eq. 0 cannot be gapped by fermion bilinear 
terms, for any flavor number n. 

We emphasize that the above Dirac fermions on the 
surface of a 3D TCI cannot be realized in any 2D sys¬ 
tem with mirror symmetry, as expected for symmetry- 
protected topological phases in general. According to the 
Hamiltonian ([Tj) , the surface states with k x = 0 within 
a given mirror subspace are chiral as they all move in 
the same direction PH- In contrast, in any 2D system 
single-particle states within a mirror subspace cannot be 
chiral (this is demonstrated with a 2D lattice model in 
Sec. [V] of the Supplemental Material [55]). 

17(1) Higgs phase and Z 2 topological order. In this 
work, we study interacting surface states of TCIs with 
n = 4. Starting from four copies of Dirac fermions in 
the noninteracting limit, we will introduce microscopic 
interactions and explicitly construct a Z 2 topologically- 
ordered phase on the TCI surface, which is gapped and 
mirror symmetric. 


Our construction is inspired by the work of Senthil and 
Fisher [26] and Senthil and Motrunich m on fractional¬ 
ized insulators. We construct on the surface of an n = 4 
TCI a Higgs phase with an xy-order parameter (b) ^ 0, 
which is odd under the mirror symmetry and gaps the 
Dirac fermions. Next, we couple these gapped fermions 
to additional degrees of freedom a y that are introduced to 
mimic a 17(1) gauge field. This gauge field a y plays three 
crucial roles: (i) the coupling between matter and a y re¬ 
stores the otherwise broken 17(1) symmetry and, thus, 
the mirror symmetry along with it; (ii) the Goldstone 
mode is eaten by the gauge boson and becomes massive; 
(iii) since the xy -order parameter carries 17(1) charge 2, 
the 17(1) gauge group is broken to the Z 2 subgroup in 
the Higgs phase. Because of these properties, the Higgs 
phase thus constructed is a gapped and mirror-symmetric 
phase with Z 2 topological order. 

We now elaborate on the construction (details of this 
construction can be found in Sec. |III| of the Supplemental 
Material [25]). First, we relabel the fermion flavors A = 
1,...,4 using a spin index s =t,4- and a 17(l)-charge 
index a = ± (unrelated to the electric charge). We take 
fermion interactions that are invariant under both the 
SU(2) spin rotation and the 17(1) rotation 

1/(1) : i/j as e io fy oa , a = ± (3) 

Moreover, we introduce a boson field b(x,y) that carries 
ll(l)-charge 2 and is odd under mirror symmetry, 

1/(1) : 6-)- e l29 b, M : b(x,y) —»■ —b(—x,y), (4) 

and couple this boson to the massless Dirac fermions as 
follows 

Hb f = Vbm 3 T~ b a z ip bs + H.c. (5) 

When these bosons condense, ( b ) ^ 0 spontaneously 
breaks both the 1/(1) and mirror symmetry, and gappes 
out the fermions. 

Finally, we introduce another boson vector field 
a y (x,y), which couples to b and ip as through minimal 
coupling. An effective theory of this system has the fol¬ 
lowing form, 

£ = - iipla^id^ + iapT^ips + (&'i/>]r+ <r z if> a + H.c.) 

+ - 2 *a M )&| 2 + r\ b\ 2 + u\b \ 4 + F^F^, ^ ^ 

where the matrices a 0 = 1 , a x = cr y and a v = a x . 
Furthermore, we add to the effective action an inter- 
ation term UN 2 , where N = r z ip + 2 Ub — V • E 
(Ei = F 0i = doai — diao is the electric field strength). 
In the limit of U —> oo, this enforces the local constraint 
N = 0. As a result, the bare fermion i/j s and boson b 
are no longer low-energy excitations, since adding them 
to the ground state violates the constraint N = 0 and 
costs an energy U . Therefore, in the low-energy effective 
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model ip s and b must be screened by the gauge field 
and become quasiparticles tp as = V’ a S e xa9 and b = be 219 , 
where the operator e m9 creates n gauge charge of and 
restores the constraint N = 0. In terms of these quasi¬ 
particles, the effective theory becomes 

£ = - + ia^r^s + a z ^ s + H.c.) 

1 _ (7) 

+ —\(d M - 2ia y )b\ 2 + r\b\ 2 + u\b\ 4 + F^F^. 

^9 

Furthermore, a U{ 1) gauge symmetry emerges in the low- 
energy Hilbert space defined by the local constraint N = 
0 [28]. Specifically, the constraint is the Gauss law and it 
restricts the low-energy Hilbert space to states that are 
invariant under the gauge transformation 

U# : -»■ e^ips, b^be 2^c|, , ->• - d y (j). (8) 

In this effective theory with the emergent 17(1) gauge 
field, condensing the boson b no longer breaks the global 
U( 1) and the mirror symmetries, as it instead breaks the 
U( 1) gauge symmetry to Z 2 . Naively, the mirror sym¬ 
metry maps (b) to —(b). However these two symmetry 
breaking vacuua are equivalent because they are related 
by the gauge symmetry transformation U n /2 '■ b ► —b. 
This restoration of mirror symmetry becomes clearly 
manifested if we assume that b and tp s transform projec- 
tively under mirror symmetry with the additional U(l) 
gauge transformation U n / 2 , 


M : 4> s (r) -)• it z <g> cr x $ s (r'); b(r) ->■ b(r'). (9) 


This Higgs phase obtained by condensing charge-2 b 
field indeed has a Z 2 topological order when the number 
of Dirac fermions is n = 4 mm- This can be under¬ 
stood by identifying the Bogoliubov quasiparticle ip and 
vortices as the anyons e, m, and e (see Sec. [IT] of the 
Supplemental Material [25] for the definition of the nota¬ 
tion) in the Z 2 topological order. ip becomes the e anyon 
as both are fermions. As the Higgs field gaps out four 
Dirac fermions, there are four Majorana fermions, or two 
complex fermion zero modes, in each vortex core. Hence, 
there are two types of vortices whose core has even or 
odd fermion parity, respectively. In the case of n = 4, it 
can be shown that the vortices carry Bose statistics (see 
Sec. IV of the Supplemental Material [25] for details), 
and they are mapped to the m and e anyons in the Z 2 
topological order, respectively. 

Mirror-symmetry fractionalization. Now we consider 
how the mirror symmetry acts in the Z 2 spin liquid phase 
described by Eq. ([7|. In this effective theory, the ip field 
is the fermionic anyon e. Equation ([ 9 ]) implies that it 
carries M 2 = — 1. 


Next, we consider how the mirror symmetry acts on 
the m anyon, which is a vortex of the Higgs field where 
all core states are empty. Since the mirror symmetry pre¬ 
serves the Higgs field (b) but maps x to —x, it maps a 
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FIG. 1. (a) A vortex and an antivortex. The direction of 

the arrow represents the phase of the Higgs field (b). The 
dotted line is the mirror axis, (b) Illustration of the fermion 
spectrum flow from the left to the right as we create a vortex- 
antivortex pair from the vacuum and move them far apart. 
In this process, four vortex core states are separated from the 
bulk spectrum, two from the conducting band and two from 
the valence band, and become degenerate zero modes. As 
illustrated in the inset, the core states are twofold degenerate 
with spin s = ±1. 


vortex to an antivortex. Therefore we consider a mirror- 
symmetric configuration with one vortex and one antivor¬ 
tex, as shown in Fig. [lja). The symmetry fractionaliza¬ 
tion of M 2 = ±1 can be detected from the M parity of 
the fermion wave function with such a vortice-antivortex 
pair |30l : for two bosonic vortices, the mirror parity is 
equal to M 2 . 

From Eq. 0, we get the fermion Hamiltonian 

H = vtpl ( k x o y - k y a x ) ^ + ^((6)r + + (b)*T~)cr z 'ii> s . 

(10) 

It has a particle-hole symmetry S : if) —> a z il> which maps 
H to SU5 = —H. This implies that its spectrum is sym¬ 
metric with respect to zero. Assume that the dimension 
of the whole Hilbert space is 4IV; there are 2 N states 
with positive energy and 2 N states with negative energy. 
For the vortex configuration in Fig. [l|a), there are four 
complex zero modes, two from each vortex core, which 
are all unoccupied. Therefore, excluding these four states 
there are 2 N — 2 states with negative energy, which are 
all occupied in the fermion wave function. 

Next, we consider the mirror eigenvalues of these 
2N—2 occupied fermion states. Since ip carries M 2 = — 1, 
each state has mirror eigenvalue Xm = ±i. Because both 
H, S and M are diagonal in pseudospin s =f, all oc¬ 
cupied states are pseudospin doublets, and two states 
in each doublet have the same Am- Hence, the mirror 
eigenvalue of all occupied states, organized as IV — 1 dou¬ 
blets, is (—l)^ 1 = —1. Therefore, the wave function of 
two empty vortices is odd under mirror symmetry, which 
implies that the e particle has the symmetry fractional¬ 
ization M 2 = —1. Combining the results that both e 
and e carry M 2 = —1, we conclude that the m anyon 
also has M 2 = —1 (see the discussion in Sec. [h] of the 
Supplemental Material [25]1. 

In summary, the gapped Z 2 surface state we con¬ 
structed has an anomalous mirror-symmetry fractional¬ 
ization that both types of anyons carry M 2 = — 1, which 
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cannot be realized in a genuine 2D system. For com¬ 
parison, it is shown in Sec. [TT] of the Supplemental Ma¬ 
terial [25] that applying the same construction to a 2D 
lattice model results in a Z 2 phase with a different mirror- 
symmetry fractionalization which is not anomalous. 

Mirror anomaly. The anomalous crystal symmetry 
fractionalization presented in the surface topological or¬ 
der implies a symmetry-protected topological degener¬ 
acy associated with the edges of the surface topological 
ordered region. This mirror anomaly is a remnant of 
the anomalous surface fermion modes in the free-fermion 
limit. To see this, we consider the setup presented in 
Fig-! in which the Z 2 surface topologically ordered state 
is terminated at two edges symmetric with respect to the 
mirror plane, by two regions with opposite ( b ) = ±1 on 
either side of the mirror plane, respectively. 



x = 0 


FIG. 2. Two mirror-symmetric edges of a Z 2 surface topolog¬ 
ical order. The mirror symmetry maps x to — x with respect 
to the mirror plane marked by the red disk at x = 0. The 
surface topological order marked by the shade terminates at 
edges against two ordered regions with (b) = ±1, respectively. 

This setup itself does not break the mirror symmetry, 
and all local excitations can be gapped everywhere on 
the surface. In particular, since the Z 2 topological order 
is not chiral, its edge can be gapped out by condensing 
either e or m anyons on the edge m- The edges next to 
an ordered phase with (b) ^ 0 are e edges, as condens¬ 
ing e breaks the global U(l) symmetry. A Z 2 spin liquid 
state on an infinite cylinder has four degenerate ground 
states |'F a ), each has one type of anyon flux a = 1, e, to, e 
going through the cylinder. On a finite cylinder with two 
e edges, only I'Ll) and |'k e ) remain degenerate, because 
adding an to or e anyon on the edge costs a finite energy. 
In a generic Z 2 state, this degeneracy can be further 
lifted by tunneling an e anyon between the two edges, 
H t = Ae^e^ + H.c., where e' L R creates two e anyons on 
the two edges, respectively. However, the e anyon carries 
M 2 = —1; therefore, the tunneling term H t is odd un¬ 
der mirror and, thus, forbidden by M. As a result, this 
twofold topological degeneracy is protected by the mir¬ 
ror symmetry even in the limit of L —> 0. This argument 
is formulated using the effective edge Lagrangian in the 
Supplemental Material [ 25 ] ■ 

In the limit of L —> 0, this topological degeneracy be¬ 
comes a local degeneracy protected by the mirror sym¬ 
metry. Therefore if the Z 2 topological order is killed by 
collapsing two gapped edges, the ground state is either 
gapless or mirror-symmetry breaking, and this cannot 
be avoided regardless of edge types because all types of 
anyons have M 2 = — 1. This topological degeneracy re¬ 


veals the anomalous nature of this mirror-symmetry frac¬ 
tionalization. Furthermore, if we collapse two gapless 
edges of the Z 2 state, the edges remain gapless because 
the anyon tunneling is forbidden by M. Hence we get a 
gapless domain wall with central charge c = 1 + 1 = 2, 
which recovers the edge with four chiral fermion modes 
in the aforementioned free-fermion limit. This is ex¬ 
plained in more detail in Sec. |VI| of the Supplemental 
Material [25] . 

Conclusion. In this Letter, we show that the surface of 
a 3D mirror TCI with mirror Chern number n = 4, con¬ 
taining four gapless Dirac fermion modes in the free limit, 
can be gapped out without breaking the mirror symme¬ 
try by a Z 2 topological order. This surface Z 2 topological 
order has an anomalous mirror-symmetry fractionaliza¬ 
tion in which all three types of anyons carry fractional¬ 
ized mirror-symmetry quantum number M 2 = —1, and 
such a topological order cannot be realized in a purely 
2D system. 

Our finding also puts constraints on possible ways to 
fractionalize the mirror symmetry in a 2D Z 2 quantum 
spin liquid [23] [321. The result of this work indicates 
that the combination that both the e and to carry the 
fractionalized M 2 = — 1 is anomalous and cannot be re¬ 
alized in a 2D Z 2 spin liquid. Furthermore, our result 
can be easily generalized to also rule out the combina¬ 
tion that the e anyon carries spin-| and to anyon carries 
M 2 = — 1 |33], because if e carries M 2 = +1, we can 
define a new mirror symmetry M' = Me l7rS , for which 
both e and to carry (M') 2 = —1 and, therefore, this 
combination is also anomalous. In summary, our finding 
implies that the vison must carry M 2 = +1 in a Z 2 spin 
liquid where the spinon carries a half-integer spin. 

We thank Chen Fang and Senthil Todadri for in¬ 
valuable discussions. Y.Q. is supported by National 
Basic Research Program of China through Grant No. 
2011CBA00108 and by NSFC Grant No. 11104154. L.F. 
is supported by the DOE Office of Basic Energy Sci¬ 
ences, Division of Materials Sciences and Engineering, 
under Award No. DE-SC0010526. This research was 
supported in part by Perimeter Institute for Theoretical 
Physics. Research at Perimeter Institute is supported 
by the Government of Canada through Industry Canada 
and by the Province of Ontario through the Ministry of 
Research and Innovation. 


[1] M. Hasan and C. Kane, Rev. Mod. Phys. 82, 3045 (2010) 

[2] X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 

( 2011 ) 

[3] J. E. Moore, Nature 464, 194 (2010) 

[4] A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. 
Ludwig, Phys. Rev. B 78, 195125 (2008) 

[5] L. Fidkowski and A. Kitaev, Phys. Rev. B 83, 075103 

( 2011 ) 





5 


[6] S. Ryu and S.-C. Zhang, Phys. Rev. B 85, 245132 (2012) 

[7] H. Yao and S. Ryu, Phys. Rev. B 88, 064507 (2013) 

[8] X.-L. Qi, New J. Phys. 15, 065002 (2013) 

[9] L. Fidkowski, X. Chen, and A. Vishwanath, Phys. Rev. 
X 3, 041016 (2013) 

[10] C. Wang and T. Senthil, Phys. Rev. B 89, 195124 (2014) 

[11] L. Fu, Phys. Rev. Lett. 106, 106802 (2011) 

[12] Y. Ando and L. Fu, Annu. Rev. Condens. Matter Phys. 
6, 361 (2015) 

[13] A. Alexandradinata and B. A. Bernevig, 
irXiv:1411.3338 

[14] T. H. Hsieh, H. Lin, J. Liu, W. Duan, A. Bansil, and 
L. Fu, Nat. Commun. 3, 982 (2012) 

[15] Y. Tanaka, Z. Ren, T. Sato, K. Nakayama, S. Souma, 
T. Takahashi, K. Segawa, and Y. Ando, Nat. Phys. 8, 
800 (2012) 

[16] P. Dziawa, B. J. Kowalski, K. Dybko, R. Buczko, 
A. Szczerbakow, M. Szot, E. usakowska, T. Balasubra- 
manian, B. M. Wojek, M. H. Berntsen, O. Tjernberg, 
and T. Story, Nat. Mater. 11, 1023 (2012) 

[17] S.-Y. Xu, C. Liu, N. Alidoust, M. Neupane, D. Qian, 
I. Belopolski, J. D. Denlinger, Y. J. Wang, H. Lin, L. A. 
Wray, G. Landolt, B. Slomski, J. H. Dil, A. Marcinkova, 

E. Morosan, Q. Gibson, R. Sankar, F. C. Chou, R. J. 
Cava, A. Bansil, and M. Z. Hasan, Nat. Commun. 3, 
1192 (2012) 

[18] Y. Okada, M. Serbyn, H. Lin, D. Walkup, W. Zhou, 
C. Dhital, M. Neupane, S. Xu, Y. J. Wang, R. Sankar, 

F. Chou, A. Bansil, M. Z. Hasan, S. D. Wilson, L. Fu, 
and V. Madhavan, Science 341, 1496 (2013) 

[19] I. Zeljkovic, Y. Okada, M. Serbyn, R. Sankar, D. Walkup, 
W. Zhou, J. Liu, G. Chang, Y. J. Wang, M. Z. Hasan, 
F. Chou, H. Lin, A. Bansil, L. Fu, and V. Madhavan, 
Nat. Mater. 14, 318 (2015) 

[20] H. Isobe and L. Fu, Phys. Rev. B 92, 081304(R) (2015) 

[21] J. C. Y. Teo, L. Fu, and C. L. Kane, Phys. Rev. B 78, 
045426 (2008). 

[22] X.-G. Wen, Phys. Rev. B 44, 2664 (1991), Int. J. Mod. 
Phys. B 05, 1641 (1991) 

[23] X.-G. Wen, Phys. Rev. B 65, 165113 (2002) 

[24] A. M. Essin and M. Hermele, Phys. Rev. B 87, 104406 
(2013) 

[25] See the Supplemental Material, which includes Refs. |34j - 
137] , for reviews of 3D TCIs, 2D Z 2 spin liquids, and their 
edge theory. 

[26] T. Senthil and M. P. A. Fisher 
(2006) 

[27] T. Senthil and O. Motrunich, 

( 2002 ) 

[28] Here the constraint N = 0 can be enforced locally be¬ 
cause the U(l) symmetry defined in Eqs. and |4| is 
an on-site symmetry operation. Particularly the U(l) ro¬ 
tation of the Dirac fermion ip is a local rotation between 
four flavors of Dirac fermions, which arise from the four 
flavors of Dirac fermions in the bulk of the 3D TCI, as 
explained in Sec. [I] of the Supplemental Material [2'5j . 

[29] M. A. Metlitski, L. Fidkowski, X. Chen, and A. Vish¬ 
wanath, arXiv: 1406.3032 

[30] Y. Qi and L. Fu, Phys. Rev. B 91, 100401(R) (2015) 

[31] M. Barkeshli, E. Berg, and S. Kivelson, Science 346, 722 
(2014) 

[32] M. Barkeshli, P. Bonderson, M. Cheng, and Z. Wang, 
rrXiv:1410.4540 

[33] M. Hermele and X. Chen, arXiv: 1508.00573 


, Phys. Rev. B 74, 064405 
Phys. Rev. B 66, 205104 


[34] Y.-M. Lu, G. Y. Cho, and A. Vishwanath, 
irXiv: 1403.0575 

[35] A. Y. Kitaev, Ann. Phys. (N. Y). 303, 2 (2003) 

[36] S. Sachdev, Phys. Rev. B 45, 12377 (1992) 

[37] A. Tanaka and X. Hu, Phys. Rev. Lett. 95, 036402 
(2005) 





Supplemental Material 


Yang Qi 1,2 and Liang Fn 3 

1 Institute for Advanced Study, Tsinghua University, Beijing 100084, China 
2 Perimeter Institute for Theoretical Physics, Waterloo, ON N2L 2Y5, Canada 
3 Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 


I. MIRROR TOPOLOGICAL CRYSTALLINE INSULATOR 

In this section we briefly review non-interacting TCIs protected by the mirror symmetry 
x —» — x. An example of such a TCI with mirror Chern number n can be constructed from 
the following low-energy Dirac fermion Hamiltonian in three dimensions 

n 3 

H = I k i T i + mr 0 )cA (1) 

A= 1 i= 1 

where the 4x4 matrices r, with i — 0,1, 2, 3,4 form a set of anticommuting Dirac Gamma 

matrices, and a is the flavor index. The fermion field ca transforms as follows under mirror: 

M : c A (x,y,z) ->• T 14 c A (-x,y, z), r 14 = ?T 1 r 4 . (2) 

It is straightforward to check that 

H = MHM~\ (3) 

i.e., the Hamiltonian (1) is invariant under mirror. As we will explain in Section II, for the 
purpose of this work, it is conceptually convenient to regard the four-components of the 
Dirac fermion in (1) as unrelated to electron’s spin. Relatedly, the mirror operation defined 
here is a Z 2 symmetry with the property M 2 = 1. In our previous works on spin-orbit 
coupled systems, the mirror operation M' acts on electron’s spin in addition to its spatial 
coordinate, which leads to M' 2 = —1. Nonetheless, one can redefine the mirror operation 
by combining M' with the U (1) symmetry of charge conservation c —>■ ic, which restores the 
property M 2 = 1. 

As discussed in the main text, the mirror Chern number n is defined as the difference 
n = n. |_ — n_, where n± are the the Chern numbers of single-particle states with mirror 
eigenvalue M = ±1 on the mirror-symmetric plane k x = 0 in the 3D Brillouin zone. Starting 
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from a trivial insulator with n + = n_ = 0, a TCI with a nonzero mirror Chern number n 
is obtained by going through a band inversion transition that corresponds to reversing the 
sign of the Dirac mass m for all n flavors in the Hamiltonian (1). The total Chern number 
is unchanged during this transition. 

II. MIRROR SYMMETRY FRACTIONALIZATION IN Z 2 SPIN LIQUID STATE 

In this section we discuss the possible ways of mirror symmetry fractionalization in a two- 
dimensional Z 2 spin liquid state. We first discuss this in the context of a bosonic system, 
then generalize the results to systems with fermions. 

The Z 2 topological order contains four types of anyons: 1, e, m, and e, denoting the 
trivial bosonic particle, the bosonic spinon, the vison, and the fermion spinon, respectively. 
Except the trivial anyon 1, the other three anyons can carry fractionalized mirror symmetry 
representation M 2 = ±1. Because the e anyon is a bound state of e and m anyons, its mirror 
symmetry fractionalization can be derived from the ones of e and m anyons as 

Ml = (4) 

where the extra minus sign appears because the e and m anyons have a nontrivial mutual 
statistics, which gives rise to the Fermi statistics of the e anyon [1, 2], Using this relation 
we can enumerate the four possible ways of mirror symmetry fractionalization in a Z 2 state, 
as listed in Table I. 


TABLE I. Four different ways of mirror symmetry fractionalization. 


No. 

e 

m 

e 

Realization 

1 

+1 

+1 

-1 

Standard toric code 

2 

-1 

+1 

+1 

Some spin liquids 

3 

+1 

-1 

+1 

Same as above 

4 

-1 

-1 

-1 

Surface of 3D TCI 


Among the four states listed in Table I, three can be realized in a purely two-dimensional 
system. In the first state the e and m anyons transform trivially under mirror symmetry, 
and this is realized in the standard Kitaev toric code model [3]. In the second state, the e 
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anyon has M 2 = —1, and this is realized in many mean field spin liquid states, for example 
the Q i = Q -2 state on the kagome lattice [4], Without additional global symmetries the e 
and m anyons are interchangable in a Z 2 topological order, so the third state is equivalent 
to the second state. 

Since the first three states can all be realized in purely 2D systems, this leaves the fourth 
state as the only possibility of a Z 2 state with an anomalous mirror symmetry fractionaliza- 
tion. Indeed, as we see in the main text, this state can only be realized on the surface of a 
3D nontrivial TCI. 

In the above discussion we considered a purely bosonic system where all physical degrees 
of freedom are bosons. However, the system we study in the main context is an interacting 
fermionic system. In this case, besides the aforementioned anyons {l,e,m, e} the system 
also has local fermion /, which carries no gauge charge or flux but obeys Fermi statistics, 
and it can form bound states with other types of anyons. Particularly the bound state ex/ 
is a boson and has a —1 mutual statistics with either e or m anyon. Therefore it is similar 
to the e or m anyon, and the set of anyons {1, e, e x /, m x /} have the same self and mutual 
statistics as the original set of {l,e,m,e}. Hence if we neglect the fermion parity, the three 
anyons e, m and ex/ becomes indistinguishable. 

The mirror symmetry fractionalization of ex f can be derived from the fractionalization 
of e anyon and the mirror symmetry representation of / fermion, 

M, 2 x/ = MlM) = -MlMlM). (5) 

Hence for the anomalous mirror symmetry fractionalization we are considering (State 4 in 
Table I), M 2 x j = — M 2 . If the mirror symmetry acts as M 2 = —1 on / fermion, the bound 
state ex f transforms trivially under M and can be condensed without breaking M. Such an 
anyon condensation confines all anyons and drives the surface to a trivial state. Therefore the 
surface topological order is only anomalous if fermion carries M 2 = +1, which is consistent 
with the fact that M 2 = — 1 does not protect topological crystalline superconductors in 3D 
(for TCI the mirror symmetry can be redefined to have M 2 = +1 because of the charge 
U(l) symmetry, as explained in the main text.) This is the reason why we choose to define 
M 2 = +1 for the fermions in this work. Furthermore, in the case of M 2 = +1, in the 
anomalous case all the three bosonic anyons e, m and ex f has M 2 = —1. Therefore we do 
not need to examine the fermion parity of the anyon excitations when studying their mirror 
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symmetry fractionalization. 

III. CONSTRUCTION OF A SURFACE Z 2 TOPOLOGICAL ORDER. 

In this section, we review the detailed steps of obtaining a Z 2 spin liquid state in Eq. (7) 
of the main text from the surface effective theory in Eq. (6) of the main text. 

The effective theory in Eq. (6), containing the surface Dirac fermions i[) as , the scalar 
boson held b and the U(l) vector held o /t , only has a global U(l) symmetry. Specifically, 
the effective Lagrangian is invariant under the U(l) symmetry transformation where is 
invariant the ^ as and b transforms as in Eqs. (3) and (4) of the main text, respectively. We 
emphasize that in this U(l) symmetry transformation the phase factor e 1 '^ is a space-time- 
independent constant. 

Next, we promote this U(l) global symmetry to a U(l) gauge symmetry by introducing 
an interaction term UN' 2 , 

C = -i^la^(d^a^T z )^ s +(b^lT + a z ^ s +hx)+^\(di l -2iai l )b\ 2 +r\b\ 2 +u\b\ A -\-F IJlv F 111 ' +U N 2 , 

^9 

( 6 ) 

where N is the U(l) charge density N = + 2 b'b -X7-E. Here E is the held strength 

of the vector held Ei = F 0i = doai — d l ao- In the limit of strong interaction: U —* +oo, 
the last term in Eq. (6) dominates. At the leading order of the perturbation theory, we 
only consider the U term, and the system has degenerate ground states described by held 
conhgurations satisfying the constraint N = 0 everywhere. Higher orders of perturbation 
theory can then be described as an effective held theory with the constraint N = 0. 

The local constraint N = 0 can be understood as an emergent gauge symmetry. First, 
this constraint is the same as the Gauss law, 

V • E = + 2 b ] b. (7) 

Second, under the U(l) gauge transformation 

17(1)0 : ip s -A b^be 2l4> , -A - <9 M 0, (8) 

the phase of the local Hilbert space changes by e lN ^. Therefore, the low-energy Hilbert 
space satisfying the constraint of N = 0 consists of conhgurations that are invariant under 
the U(l) gauge transformation in Eq. (8). In other words, enforcing the constraint of N = 0 
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in the large-?/ limit restricts the low-energy dynamics of the model to an effective Hilbert 
space with an emergent U(l) gauge symmetry. 

Another way of understanding this emergent U(l) gauge symmetry is to notice that 
the bare excitations ij}\ and b' violates the local constraint N = 0, and creates excitations 
with U(l) charge one and two, respectively. Because of the U- term, these excitations carry 
interaction energy U and 4 U, respectively, and are therefore excluded from the low-energy 
Hilbert space in the large-?/ limit. In other words, a fermion excitation must be screened 
by gauge held fluctuations and form a quasiparticle excitation ipl = that satisfies the 

constraint of N = 0, where the operator e ld creates a unit of gauge charge: [e ld ^ x \ V -E(x')\ = 
6 2 (x — x') and restores the constraint in Eq. (7). Similarly a boson excitation ?/ must be 
screened by e 2ie and become a quasiparticle excitation ?/ = b^e 2ld . These fermionic and 
bosonic quasiparticle excitations create excited states within the low-energy Hilbert space 
defined by the constraint of N = 0, and obey the gauge transformation in Eq. (8) and Eq. (8) 
of the main text. They are used to construct the low-energy effective theory in Eq. (7) of 
the main text. 

Starting from the low-energy effective theory with the emergent U(l) gauge symmetry, 
condensing the bosonic quasiparticle b gives rise to a Higgs phase with no mirror-symmetry- 
breaking, no gapless excitations and a Z 2 topological order, in contrast to a different Higgs 
phase obtained by condensing the bare boson b in the original free-fermion model, where 
the mirror symmetry is broken. In the main text, it is explained that the b held is odd 
under mirror reflection [see Eq. (4) of the main text]. Hence, a condensate with ( b} ^ 0 is 
mapped to a different condensate with —(b) under the mirror symmetry, indicating that the 
Higgs phase of (b) ^ 0 breaks the mirror symmetry. However, the condensate of (b) does 
not break the mirror symmetry. Although naively mirror symmetry maps b to —b, these 
two held configurations are related by a gauge transformation of U n / 2 that maps —b back 
to b, and they represent the same physical state in the low-energy Hilbert space defined 
by the constraint iV = 0, as a result of the emergent U(l) gauge symmetry. Therefore the 
condensate of ( b) ^ 0 nolonger breaks the mirror symmetry. This is more clearly manifested 
by the redefinition of the mirror symmetry action in Eq. (9) of the main text. 

Moreover, the condensation of b gaps out all excitations and creates a Z 2 topological 
order. On one hand, the condensation generates a mass gap for the fermionic quasiparticlc 
ijj through the Yukawa coupling (bxjjlr + a z 'ip s + h.c.'). On the other hand, it also generates a 



6 


mass gap for the gauge boson through the Higgs mechanism. Since the Higgs held b carries 
two units of the gauge charge [see Eq. (8) of the main text], it breaks the emergent U(l) 
gauge symmetry to a Z 2 gauge symmetry. This Z 2 gauge symmetry is generated by the 
remaining gauge transformation U n , which maps ip s to —ip s . Therefore, in the Higgs phase, 
the fermionic quasiparticle ip s carries a unit of Z 2 gauge charge, and it couples to a Z 2 gauge 
held, which is the reminant of the emergent U(l) gauge held. As shown in Sec. IV, this 
Higgs phase has the Z 2 topological order. 


IV. TOPOLOGICAL ORDER IN THE HIGGS PHASE 

In this section, we show that the Higgs phase obtained in the main text after condensing 
(b) has the Z 2 topological order. Naively the condensation of the Higgs held breaks the U(l) 
gauge symmetry to Z 2 and therefore the remaining gauge fluctuation is described by a Z 2 
gauge theory. However, integrating out the gapped Dirac fermions can produce nontrivial 
Berry phases that alters the statistics of the vison excitations in the Z 2 gauge theory. Here 
we show that in our case the visons still have trivial bosonic statistics and the Higgs phase 
has the Z 2 topological order. 

After the Higgs condensation the Bogoliubov quasiparticles and the vortices become 
anyon excitations. First, the Bogoliubov quasiparticle ip is a fermion and carries a Z 2 
charge. Therefore it can be identified as the e anyon. Next, we consider the vortices of 
the b condensate in this Higgs phase. We note that such vortices host four Majorana zero- 
energy core states in the vortex core, one from each flavour of Dirac fermion, and these four 
Majorana zero modes can be grouped into two complex fermion zero modes. Hence there 
are four degenerate core states, labeled by the occupation of these two complex zero modes 
I 711712 ). This four-fold degeneracy can be lifted by fermion interaction, after which the two 
states with even fermion parity (|00) and 111)) and the two states with odd fermion parity 
(|01) and 101)) are still degenerate. We first consider the brading statistics of the empty 
vortex 100). To do that we define an antiunitary transformation 

T : ip s —>■ ir z 0 icr y ip s ; b —> b, (9) 

which can be viewed as an analog of the time-reversal symmetry. It is easy to check that the 
effective Lagrangian of the Dirac fermions with s =t and s =j. in Eq. (7) is invariant under 
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T separately. This implies that their contribution to the Berry phase of exchanging two 
empty vortices must be ±1. Therefore combining the contributions from s =)\|, the total 
exchange phase is (±1) 2 = +1, indicating that the empty vortex is a boson. Moreover, it 
carries a n flux of the emergent U(l) gauge field (because b carries charge-2), and therefore 
has a —1 mutual statistics with the ip fermion. Hence the empty vortex can be identified 
as the e anyon. Furthermore, the half-filled vortex with one fermion mode occupied can 
be viewed as the bound state of an empty vortex and a Bogolinbov quasiparticle ip, and 
therefore corresponds to the m anyon. In conclusion, in the case of n — 4 the Higgs phase 
has the Z 2 topological order, and this is in contrast to other more exotic topological orders 
obtained in the case when n is not a multiple of four [5, 6]. 


V. MIRROR SYMMETRY IN THE 2D SQUARE LATTICE vr-FLUX MODEL. 

If mirror symmetry is not broken, a 3D TCI discussed in this work has an anomalous sur¬ 
face that either has gapless Dirac fermions or has a topological order in which anyons carry 
anomalous mirror symmetry fractionalization. It is worthwhile to compare these anomalous 
TCI surface states with the 7r-flux model on the square lattice. This genuine 2D model 
has Dirac fermions that can be gapped out without breaking mirror symmetry or adding 
interactions. Furthermore, if a Z 2 topologically ordered state is constructed using the same 
approach as described in the main text, the resulting state has mirror symmetry fractional¬ 
ization that is not anomalous. 

The 7r-flux model on a 2D square lattice is known to host 2D Dirac fermions at half¬ 
filling. By linearizing the tight-binding model and keeping track of the mirror symmetry, 
we find [7, 8] the following Dirac Hamiltonian for the 7r-flux model 

H 2 d = v ^ ^2 1 Ph( k x<7v ~ (-1 ) a k y a x )ip as . ( 10 ) 

o=l,2 s =t,t 

where the combination of a = 1,2 and o z = ±1 denotes four lattice sites in a unit cell, and 
s labels electron’s spin [9]. Each fermion field ip as transforms under the mirror symmetry 
of the square lattice as described by Eq. (2). Comparing Eq. (10) to Eq. (1) in the main 
text, we find that unlike the case of TCI surface states, states on the 2D square lattice are 
non-chiral within each mirror sector, as can be seen from the sign factor (—1)“. In this case, 



a coupling between the two Dirac points 


i’hixi’ 2 , + h.c. 


(ii) 


opens np a gap, while preserving the mirror symmetry. 

This construction of Z 2 spin liquid states also applies to the 2D Dirac fermion effective 
model in Eq. (10) derived from the 7r-flux model. However the mirror symmetry fraction¬ 
alization in the constructed Z 2 spin is different. Because the Dirac fermions can now be 
gapped out by the coupling in Eq. (11), we modify the way the b boson couples to the Dirac 
fermion in Eq. (5) to the following form, 

H' bf = Vb^l s T~ b a x ^ bs + h.c.. (12) 

Hence in the Z 2 spin liquid state the anyons b and 0 transforms trivially without an extra 
U(l) gauge transformation under mirror symmetry, and this implies that the 0 (or the e) 
anyon carries trivial mirror symmetry fractionalization M 2 = +1. Therefore only one of the 
e and m anyons carries M 2 = — 1 and the other one still carries M 2 = +1. This crystal 
symmetry fractionalization is allowed in a purely two-dimensional system, as discussed in 
Sec. 1. 


VI. EDGE LUTTINGER THEORY 

In this section we derive the conclusions on edges of a Z 2 spin liquid state using the 
effective field theory of the edges, which has the form of a Luttinger theory. 

The edge of a Z 2 spin liquid state is described by the following single-channel Luttinger 
Liquid effective action 

1 7/ 

£edge = -d t (j)dyd - — [ K(dyO ) 2 + K~ 1 (dy4>) 2 ] ■ (13) 

where 9 and 0 describes edge fluctuations of e anyon and m anyon flux, respectively, and 
they obey the commutation relation that [<f>(y),9(y , )\ = i|sgn(?/ — y'). In a Z 2 state both 
two e or two m particles fuse into a trivial excitation and therefore the perturbations e ±2ld 
and e ±2t ^ are allowed on the edge, and in general the gap will be gapped out by either one of 
the two terms, unless the gaplessness is protected by some other global symmetries. (In our 
setup, each edge by itself is not symmetric under the mirror symmetry but mapped to the 
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other edge instead; hence the mirror symmetry puts no constraints on these terms.) Here 
in our setup, the Z 2 phase is terminated by the ordered phase ( b ) ^ 0, where b ~ e 2,e . This 
corresponds to an edge that is gapped by the perturbation e ±2ie , which is described by the 
following action, 

1 01 

Ce dge = -d t( pd y e - — [K{d y ef + K-\d y <P) 2 ] + A cos 29. (14) 

7T Z7T 

In the limit that A —> 00 , the 9 field condenses at the minima 6 = 0,7T, which is consistent 
with the order state ( b ) = +1. Similarly for the right edge, the potential term is —A cos 2d, 
giving rise to minima at 0 — f, as the ordered state on the right has (b) = —1. Putting 
the left and the right edges together, the effective Lagrangian is the following, 

c = V ( -a,o„i:i,(i, - 4 [K( 0 ,e a ) 2 + ) 

AW 21 > (15) 

+ A(cos2d L — cos2dR). 

The above Hamiltonian is invariant under the mirror transformation if we assume that 6 
field transforms under mirror as 

M : 6 l r —> 0 R L + —. (16) 

Since the e anyon is identified as e ~ e , this is consistent with the mirror symmetry 
fractionalization M 2 e = —e. 

In the effective action of Eq. (15) each 9 a has two degenerate minima, so there are four 
degenerate minimal configurations labeled by \6l9r), where 9l = 0, n and 9 R = |, 
However the states \9 L 9 R ) and \9 L + tt, 9 r + tt) are related by the Z 2 gauge transformation 
and therefore represent the same physical state. Therefore there is a two-fold ground state 
degeneracy labeled by A 9 = 9 R — 9l = ±f. 

When we bring the two edges together, in a normal Z 2 state there will be tunneling terms 
that couple the two edges together and lift the two-fold degeneracy by locking the relative 
phase of A 9. To the lowest order such tunneling term has the following form, 

C t = T ie ieR - ieL + T 2 e idR+i6L + c. c.. (17) 

Using the mirror symmetry transformation in Eq. (16) we see that the mirror invariance 
requires that T\ = rf and r 2 = 0, so the tunneling term takes the form C t = TiCOs(A0), 
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which does not lift the degeneracy between A 6 = ±|. Therefore this two-fold ground 
state degeneracy remain exact as we collapse the two gapped edges and eliminate the Z 2 
topological order, and this can only be realized at the surface of a nontrivial 3D TCI. 

Furthermore, if we collapse two gapless edges we can recover the gapless Dirac fermion 
modes in the free fermion limit. To study gapless edges, we can fine tune the parameters 
in the effective action of Eq. (15) such that the edges are gapless. Without the mirror 
symmetry, the two gapless edges will be gapped by the tunneling term T 2 e ieR+ldL in Eq. (17) 
when we collapse the two edges. However, the mirror symmetry enforces that r 2 = 0, which 
implies that the two edges remain gapless (the T\ term becomes a constant in this limit since 
Or > 0l as L —» 0). Therefore when we collapse the two edges, the Z 2 spin liquid region 
becomes a domain wall separating two gapped regions, and this domain wall hosts gapless 
modes with central charge c = 1 + 1 = 2, where the left and right edges each contribute 
c — 1. This is consistent with the conclusion that in the free fermion limit such a domain 
wall hosts four Dirac fermion modes [10]. 
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